interfiles <- 'formatted_data/stream_co2'
year <- 2022

co2_2022 <- opn_concat_all_2022(interfiles, year) 

co2_2022$co2_ppm <- (co2_2022$Voltage..V./3)*10000
co2_2022 <- co2_2022 %>%
  select(., year, shed, site, datetime, co2_ppm)

year <- 2021

co2_2021 <-opn_concat_2021(interfiles, year) %>%
  select(., year, shed, site, formatted_datetime, co2_ppm) %>%
  rename(., datetime = formatted_datetime)

all_co2 <- rbind(co2_2021, co2_2022)
co2_check_plot <- ggplot(all_co2[all_co2$shed == 'dh' & all_co2$site =='dh1' ,], aes(datetime,co2_ppm, color = year)) +
  geom_point() +
  facet_wrap (~site)

ggplotly(co2_check_plot)